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The dynamical symmetry breaking in a quasi-(l+l)-dimensional relativistic model is investigated. 
- - - The motions of particles in intrachain are described as a relativistic electron-hole gas, while the 

interchain hopping term is introduced as a Oth-component of vector in (l-l-l)-dimensions, a kind 
of chemical potential of the system. The gauge symmetry of the model is chosen as (7(1) suitable 
for a possible situation of a real substance in condensed matter physics. We consider the BCS-type 
contact interactions for the s-wave fermion-pair condensates, while employ the nonlocal interactions 
O ■ of the generalized BCS framework to generate the p-, d- and /-wave condensations in the system. 

Especially we examine the dynamical generation of a Dirac mass term and superconductivity in the 
model. The phenomenon is interpreted as metal-insulator/metal-superconductor phase transitions. 

PACS numbers: ll.lO.-z, 70.20.fg, 74.20.Rp, 74.70.Kn 
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Recently, many papers have considered the competitions/coexistences of several orders coming from particle-hole 
and particle-particle type fermion pair condensations [1-11]- For example, in quasi- (l-fl)-dimensional organic super- 
'■ conductors ( the Bechgaard salts (TMTTF)2^ ) [3-6,8-11], the competition between instabilities toward taking place 
C/3 , of superconducting states ( rf-wave singlet or /-wave triplet ) and spin/charge density waves ( SDW/CDW ) was 
discussed in the model parameter space by the renormalization group analysis [9-11]. Another interesting coexistence 
of the diagonal and off-diagonal long range orders might be found in the physics of supersolid [12]. In this case, the 
possibility of the coexistence of the diagonal ( CDW ) and off-diagonal ( superfluidity ) long range orders are discussed. 
These orders seem to have a universal relation under variations with respect to external and/or model parameters 
^ > like chemical potential, temperature, coupling constant for interaction between particles: Various substances shows 
Q ' superconductivity, and in their phase diagrams, there are several orderings, CDW, SDW, Mott-Hubbard insulator, 
O . so forth, in the neighbor of the superconducting states. The competition between fermion-antifermion ( for example, 
SDW, CDW, mesons, dynamical chiral symmetry breaking ) and fermion- fermion ( diquarks, superconductivity ) 
I ] condensations are observed/considered under various situations in interacting Fermi systems, in condensed matters 
^ ■ physics, nuclear many-body systems, quark matters [13,14], so forth. This feature exists both in the description of 
^\ ' phenomenological BCS ( Bardeen, Cooper, Schrieffer ) theories [15-20] and theories on mechanisms of the origin of 
• superconductivity/superfluidity [2,8-11], because the competition of these orders has the origin in the different types 
C ' of excitations of interacting Fermi systems, particle-hole or particle-particle types. Quite a large part of physical 
' I , properties of interacting Fermi systems of condensed matters will be determined by the existence of the Fermi surface 
CNJ ' or the gap ( mass ) generation at the Fermi level. There are mainly three mechanisms of the mass generation: the 
Higgs, the BCS-NJL ( Nambu and Jona-Lasinio ) [21] and the Kaluza-Klein mechanisms. Especially in the situation 
of interacting Fermi systems of condensed matters, the BCS-NJL mechanism is the suitable method to describe the 
physical property of the systems. In this paper, we examine physical properties of a quasi-(l-|-l)-dimensional rela- 
tivistic model under competition/coexistence of several types of chiral condensations and BCS and generalized-BCS 
superconductivities [22-26] . 

This paper is organized as follows. In Sec. II, we introduce our model Lagrangian, discuss its structure and 
' symmetries. In Sec. Ill, the group theoretical consideration of the order parameters they are examined in this work 
is presented. The effective potential and the gap equations are obtained in Sec. IV. For generating not only s-wave 
condensations but also p-, d- and /- wave pairings of both fermion-antifermion and fcrmion-fermion types, we prepare 
the generalized BCS formalism. Finally in Sec. V, we give the summary of this work. The numerical part will be 
published as the part II of this study. 
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II. THE QUASI-(l+l)-DIMENSIONAL RELATIVISTIC MODEL 

We consider quasi-one-dimensional electron system [3,4]. The system has three spacetime coordinates x = 
{xq,xi,X2)- One-dimensional L chains are stacked in parallel with a fixed periodicity. The thermodynamic limit 
is taken by i ^ oo in this case. X2-coordinate is discretized, X2 = £ N, 1 < n^; < L, and denotes the site Ux- 
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Particles ( electrons and holes ) behave as a homogeneous Fermi gas toward xi ( intrachain ) direction, while they 
move to X2-direction by the interchain nearest neighbor hopping with a parameter t of the tight-binding picture [27]. 
t ( G R ) corresponds to the overlapping integrals of atomic orbitals between two nearest chains. In the vicinity 
of the Fermi energy, energy-momentum relation will be approximately described by the linear-dispersion relation 
±fF|A:i| ( vp; the (l-l-l)-dimensional Fermi velocity ). The energy spectra of particles in such a system are given as 
follows [4,9-11]: 

£^-i-(fe) = ±VF\ki\ — 2tcosk2- (1) 

Here, k = (fci, fc2). We will examine the phenomenon of dynamical symmetry breaking in the quasi-one-dimensional 
electron system. For this purpose, we employ the following Lagrangian density of a Dirac model: 

C = £em+£o+£? + £?'-'^ + £" + £s, (2) 

1 

T 

£o = i>a{xo, xi, n^) W{iDQ + vpkp -t} + ivpl^Di - u] Va(a;o, xi,n^ + 1) 



+'ipa{xo, xi,n^ + 1) \-f°{iDo + vphp -t} + ivpl'^Di - u\ e-*'=^^("''+i'"=^ Va(a;o, xi, n^), (4) 



^ = -^\{^a{x)^cc{x)f + {^cc{x)i'^r,^cc{x)f = ^{^al^Mi^ffJi^i'e), (5) 

r.f-'^ = ^ J2 J rfa;^,da;;rfj/orfyidy^dj/;V^„(a;)V'a(y)V(a;,y;a;',t/')Vi0(a;')V'/3(y')> (6) 

ny,nx>,nyi 

jCf-""" = ^ j dx'odx[dyodyidy'ody[Mx)€{y)mx,y;x',y')i^J{x')My')^ (7) 



£b = -ffjMBB-Via7°^V'/3, (8) 

where, the definitions of fields and convariant derivatives are given as follows: 

Ftiu = dij,A^-d^Aij_, tjja = {Ra{x),La{x))'^, a,/3=T,i, Do = do-ieAo, Di = di - ieAi. (9) 

£em and £o are the kinetic terms of electromagnetic and matter fields, respectively. £f is a four-fermion contact ( 
s-wave ) intrachain interaction between particles. Go is a bare coupling constant with mass dimension [mass]~^. On 
the other hand, £"'~'^ ( C^^~°'^ ) has the nonlocal interaction between particles V [W ), can include both inter- 
and intra-chain interactions. It is prepared for generating p-, d- and /- wave fermion-antifermion ( fermion-fermion 
) pair correlations, ipa is a two-component Dirac field, has the (iso)spin index a as an isodoublet, replicated into 
an TV-flavor field for taking the large- iV limit of expansion smoothly. Ra{x) ( La{x) ) describes particles move 
toward the right ( left ) direction. We should note that the Dirac field xpa is defined in (2-|-l)-dimensional spacetime, 
thus the mass dimension of the field is determined through the dimensional analysis in (2-|-l)-dimensions ( not in 
(l-l-l)-dimensions ): Va carries the mass dimension [mass]^. On the contrary, the definition of gamma matrices are 
chosen as the ordinary (l-l-l)-dimensional field theory [28,29]: 

n / 1 \ „i _ _ ^ -1 \ _ ( 1 



They satisfy the following Clifford algebra relations {7^,7'^} = '2.g^^ ( we will use the Lorentz signature g^^ = 
diag(l,— 1) throughout this paper. No euclidization will be needed ). The charge conjugation matrix C = 7^ ( 
= —7^ ) satisfies C^^C~^ = —7^^. 75 anticommutes with C in (l-l-l)-dimensions. £0 has couplings between 
matter V-'a, ''/'a and U{1) gauge fields (^ = 0,1,2,3). The local U{1) gauge invariance will be maintained under 
the following gauge transformation law [28]: 

A^{x) ^ A^{x) + U^ev{x), (m = 0,1,2,3), ^ e-'^^(^)Va(a;), 

A2{nx + l,nx)= j dsA2{s), ^2(rix + ^ ^2(nx + li'n^cc) + -( 6'y(nx + 1) - (11) 

Tlx ^ 

In the integral for the definition of j42(n2, + l, n,,,), it should be imderstood that the integration parameter ,s continuously 
runs from the coordinate of site to site + 1. The interaction between the charge density and the vector potential 
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^oV-'T^V' is already included in £o- 'Cs is the Zeeman term which would give the Pauli spin polarization, 9,7 is the 
5-factor, and jjB is the Bohr magneton. B = {Bi, B2, -B3) = —eij^djA^ ( £123 = 1, k = 1, 2, 3 ) are magnetic fields. 
The spin polarization has been introduced as a 0th component of vector in (l+l)-dimensional Lorentz symmetry. 
Hence, the Zeeman term acts as different effective chemical potentials of different spin eigenstates. In the definition 
of C, we have chosen that the third component of can couple with matter only through the Zeeman term, tpa 
and 'ipa live in (2+l)-dimensional spacetime, while the magnetic field B has three components and defined in (3+1)- 
dimensions: The gauge fields can propagate in (3+l)-dimensional spacetime. Throughout this paper, it is enough 
for us to evaluate several quantities under the (2+l)-dimensional formalism. Hereafter, in Cq are omitted, and we 
regard B as external fields ( not a dynamical degree of freedom ). The non- vanishing B fields break the time-reversal 
invariance of the system. In our system, the Fourier transform of ipa becomes 

Mxo,xun,) = J Mko,ki,k2)e-"''"'''+"'^^^+"'^^-, (12) 

The domain for the integration — tt < k2 < n corresponds to the first Brillouin zone in ^2 direction. Hence, the kinetic 
term obtains the following Fourier transformed expression: 

^0 = '>Pa[7°{ko + vpkF + IJ-{k2)} + vpj^ki - m{k2))tpa, m{k2) = '2ucosk2, ii{k2) = -2tcosk2- (14) 



Let us examine the symmetry of our model Lagrangian. Co + £f has the U{1) global gauge symmetry under the 
transformation 

V-a ^ e'^^V'a, (15) 

while it has the (l+l)-Lorcntz symmetry at Wi?fci?+/i(A'2) = 0. Of course, the 0(2, l)-Lorentz symmetry ( including the 
rotational symmetry in (a;ia;2)-plane ) in (2+l)-dimensions is explicitly broken in our model. From the viewpoint of 
(l+l)-dimensional Lorentz symmetry, we notice the fact that, we should choose the hopping term of the tight-binding 
picture as a Oth-component of vector, not a scalar, for recovering the spectra in Eq. (1) under a special condition 
of model parameters. Hence fi{k2) acts as an alternating chemical potential to the system. This fact becomes more 
clear when we diagonalize the Lagrangian to obtain the quasiparticle excitation energy spectra, which we will do in 
Sec. IV. To compare the effect and/or roles of /i(fc2) in the physical property of C we also added an alternating Dirac 
mass m{k2) as a scalar in (l+l)-dimensions into Cq. m{k2) is chosen so as to alternate in A;2-direction. In terms of 
R and L, Cq is expressed as 

Cq = tr|i?;^(«i9o + wf^i + vphp + lJ-ik2))Ra + L^iido - ivpdi +VFkF + Kk2))La - m(fc2)[-Ri^a + -^^a-Ralj- (16) 

Here, tr denotes the trace over the spin ( not spinor ) space. We have observed that, the mass term of m{k2) 
mixes modes of and ( positive and negative chirality sectors ): The global U{1) chiral symmetry under the 
transformation 

Va^e'T^^^Va, Ra^e'^^Ra, La^e-'^^La, e R, (17) 

is explicitly broken in Cq at m{k2) 7^ 0, while it is symmetric at 771(^2) = 0. The hopping term /x(A;2) is neutral under 
the chiral transformation and thus it does not break the symmetry. The Lagrangian £0 + ■Cf has the discrete chiral 
symmetry under 

tpa 75V'a, m(fc2) ^ -m{k2). (18) 

The interaction Cj has attractive channels for both particle- hole and particle-particle types at Go > 0, and this 
fact will be shown later by employing the Fierz transformation. We can consider the origin of Cj coming from 
the particle-phonon coupling as similar to the case of BCS theory [15]. A fermion-fermion repulsive interaction in 
intrachain ( with neglecting a possible long-range interaction ) similar to the case of the Hubbard interaction can 
also be treated by the contact interaction at Go < [28]. Throughout this paper, we set aside the question on 
the mechanisms of the origins of the interactions £f , Cf''^, and Cf'"^, and examine our model by the attitude of 
the BCS-NJL ( Nambu. .lona-Lasinio ) theory [15-26]. The consideration of our method is similar to the important 
works of Yamaji [5,6]. Yamaji examined the possibility of coexistence of SDW coming from the Peierls instability and 
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singlct/triplct superconductivity in a quasi- (l+l)-dimensional system, by employing an unconventional ( = generalized 
) BCS interaction into the quasi-(l-|-l)-dimensional Hubbard model. ( His conclusion is that there is no coexistence 
of SDW and singlet / triplet superconductivity. ) Our model could be classified as a relativistic version of the Yamaji 
theory, though we do not consider a possible Peierls-type instability of a mode of a finite momentum vector Q arising 
from a nesting of the Fermi surface through a response function of our model. For example, there are many papers 
to investigate the origin of effective attractive interaction in a particle-particle ( not a particle-hole ) channel for 
a superconductivity from the Hubbard-type contact repulsive interaction, through the spin fluctuation [2,8] or the 
Kohn-Luttinger mechanism [30]. In Ref. [9-11], the response functions ( namely, the mass terms of order parameters, 
a collective field is tachyonic or not ) were examined for the investigation of the instabilities toward taking place of 
ordered states. We have interest more on the physical properties of ordered states than instabilities toward ordered 
states, thus our attitude is closer to the theory of the generalized BCS theory [17-20] than theories on the mechanism 
of the origin of an attractive interaction in a superconductivity [2,8-11]. A possible anisotropy of particle interaction 
which might be arisen from the inequivalence of the direction xi and X2 ( or, ki and k2 ) is not taken into account 
in Cj at least at the classical level of our model Lagrangian. On the other hand, as mentioned above, we will use 
and C^^'^"'^ to contain possible anisotropics of interactions in our model for incorporating the generalized BCS 
framework in our theory. To discuss the symmetries of and we must specify the structure of potentials 

V and W, and this will be done in Sec. IV. 

Here, we wish to comment on the relation between our model and several (l-l-l)-dimensional field theories. If 

we omit Cem, Cf ''', C'f ^"'^ and Cb, the (l+l)-dimensional limit of our model will be obtained at t = m = 0. In 
(l-l-l)-dimensional field theories, usually the Abelian or the non-Abelian bosonization schema will be performed [31- 
37], and physical properties of the system are discussed in the context of a realization of the Tomonaga-Luttinger 
liquid. However, we do not choose this way. The (l-l-l)-dimensional relativistic model with including interaction 
terms of forward scattering is called as the Tomonaga-Luttinger model [31]. £o at f = u = corresponds to the 
kinetic term of the Tomonaga-Luttinger model. The Tomonaga-Luttinger model is exactly soluble by the method 
of bosonization [31,33], and it shows the renormalized gapless ( namely, massless ) linear dispersion E{k) = Vp\k\ ( 
Vpi A renormalized Fermi velocity ). The massless Thirring model is also solved by the Abelian bosonization, and 
it becomes a free Bose theory [32,37]. In the framework of bosonization, a large class of (l+l)-dimensional models 
become the quantum sine-Gordon theory. In that case, the effects of discrepancies of the models from the Tomonaga- 
Luttinger model have been examined by renormalization group analysis: The results in the renormalization group 
analysis show that the spectra of several models become gapless ( massless ) when the cosine term is irrelevant, 
while others have gapful ( massive ) spectra when the cosine term of the sine-Gordon theory is relevant. Co + -C/ 
at m(fc2) = l-i'{k2) = kp = coincides with the massless Thirring model, while it becomes a Hubbard-like model at 
vp = u = 0. At vp ^ t,u {t,u ^ vp ), the kinetic term ■ijj{—ivp^^di)tlj ( the hopping/mass term 'ip[—2{u+^°t) cos k2]tp 
) dominates the physical property of the motion of particles of the system. It was shown that the massive Thirring 
model of (l+l)-dimensions becomes a sine-Gordon theory, where the mass term of the Thirring model becomes the 
cosine term of the sine-Gordon theory. Physical fermion in the sine-Gordon theory is solition/antisoliton and they 
become massive under some conditions [35] . The massive Thirring model can be converted into the antiferromagnetic 
XyZ-Heisenberg chain through the Jordan- Wigner transformation [28,29,39]. The two-flavor massless Schwinger 
model is also mapped into the S = 1/2 antiferromagnetic spin chain [42]. The massless Schwinger model on the 
circle will have the chiral condensate ('0V) 7^ [42]. Moreover, the (l-l-l)-dimensional massless Gross-Neveu ( 
Thirring ) model in real space has a (5-function interaction, it was also solved by the Bethe ansatz [40]. Another 
important aspect of the (l+l)-dimensional Gross-Neveu model is that it has the asymptotic freedom [43], and it 
might be an intersting problem for quantum field theoretical point of view that how the asymptotic freedom will be 
modified in our quasi-(l-|-l)-dimensional model ( the pure-(2-|-l)-dimensional Gross-Neveu model has no asymptotic 
freedom ), though it may not be important for the low-energy phenomena of our model. The (1-f l)-dimensional 
Hubbard model describes strongly-correlated electrons where they interact by a short-range ( (5-function ) repulsive 
potential. The exact solution of Hubbard model at the half-filling case was obtained by the nested Bethe ansatz, and 
the ground state is a spin-singlet Mott insulator given by a pair condensation of particles and holes [38]. It should 
be noted that, a particle- hole transformation in the Hubbard interaction gives attractive channels between particles 
and holes. Of particular interest in this paper is the possible mechanism of the metal-insulator/metal-superconductor 
transitions by dynamically generated Dirac mass or Cooper pair condensations. We use the method of auxiliary 
fields and the large- expansion, a kind of Hartree-Fock mean field theory. We investigate the effects of the finite 
interchain hopping and the alternating mass in the dynamical generation of mass-gap and superconductivity. Due to 
the Mermin-Wagner-Coleman theorem, there is no long-range order in (1+1)- and (2+1)- dimensional field theories 
with short-range interactions [44-46]. If we evaluate a two-point correlation function, it may show the power low 
decay, and its infrared behavior is singular [29]. In the theories for examinations of phases/orders in quasi-(l+l) 
dimensional systems, usually the word "the real world is (3+l)-dimensions" has been used for the explanation on 
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this problem [5,6,9-11]. Here, we also obey this attitude. The fact proved in Ref. [35] is that, in (l+l)-dimensional 
four-fermion models, the amplitude of fcrmion-pair condensation can well develop, while the phase ( a bosonic degree 
of freedom ) of the condensation shows a power-law decay. It is just the Kosterlitz-Thouless transition [47]. Thus, at 
least on the discussion of the development of the amplitude of fermion-pair condensation, we believe we can use the 
laxge-N mean-field theory. 

We should also mention the relation between our model and other (2-|-l)-dimensional relativistic theories. The 
(2+l)-dimc>nsioiial Gross-Nc;veu model has been used for the investigations of dynamical symmetry breaking under 
various situations, under finite temperature and density, under the existence of external magnetic field, in Minkowski 
and curved spacetimes [25,48-51]. By employing the Nambu— Jona-Lasinio model and the Gross-Neveu model, the 
relativistic theory of superconductivity in (3+l)-dimensions was extensively studied until now, both for examination 
on the relativistic effects in condensed matters [22-24] and from the context of color superconductivity in QCD4 [14]. 
In Ref. [25], the relativistic method of superconductivity in (2-|-l)-dimensions was investigated by using the (2-1-1)- 
dimensional Gross-Neveu model, for the examination of superconductivity in two-band systems with honeycomb 
lattice structure, in which case the energy dispersion of electrons/holes is described by a Dirac Lagrangian. Our 
theory considered in this paper can be interpreted as a modification of the (2-|-l)-dimensional Gross-Neveu model by 
introducing the hopping terms and the nonlocal interactions. 

III. THE GROUP THEORETICAL CONSIDERATION OF THE ORDER PARAMETERS 

In this section, we consider the symmetry properties of order parameters coming from the condensations of particle- 
hole and particle-particle pairs. The main purpose of the group theoretical consideration is to specify suitable order 
parameters for our discussion of this paper. 

First, we examine the symmetry of order parameters in (l+l)-dimcnsional Lorentz ( spinor ) space. For this 
examination, the Fierz transformation of the four-fermion interaction Cj is convenient for us, and the result may also 
be helpful to handle the nonlocal interactions CJf~'^ and Cf~°'^. By utilizing the complete set {1,7'', 75} ( = 0, 1 ) 
in 2 X 2 matrix spinor space, the Fierz rearrangement in the four-fermion interaction was obtained into the following 
form: 



Here, the spin index a of V^-field is not written explicitly for simplicity. Unlike the case of (3-|-l)-dimensions, here 
we do not have any axial-vector channels like V'T'^TsV') because the matrix space is four-dimensions and the linear 
independent set for the expansion of the space will be constructed only by scalar, vector and pseudoscalar. By 
decomposing the Dirac field -0 into the left and right moving modes under our convention of the gamma matrices, 
one finds the following expressions for densities of scalar qs^ vector ( = 0, 1 ) and pseudoscalar qp ( "pion" ) of 
fermion-antifermion pair types ( "diagonal" part of the density matrix, GDW/SDW types ): 



The non-vanishing VEV ( vacuum expectation value ) of ( chiral condensation ) indicates the realization of 
spontaneous chiral symmetry breaking of the model £0 + 'Cf , while a non-zero VEV oi Qp can be called as a kind 
of "pion" condensation. It is clear from these expressions, the interaction (jp'ip)'^ includes the Umklapp processes 
B)B) LL and L^L^RR, while they will be removed by taking the combination (V'V')^ + (V'^lsV')^- The global U{1) 
chiral symmetric interaction does not have the Umklapp processes. On the other hand, fermion-fermion pair densities 
of scalar ^5, vector ( = 0, 1 ) and pseudoscalar ( "off diagonal" part of the density matrix, Gooper pair types 



= \ [(V;75C^^)(^^C-S5V') + (Vi7%5C^Vi^)(V'^C-S57V) 

-(V57'75C^f'')(V'^C-i757V) + (V^Cf^K^'^'C-V)" . 



(19) 



Qs{x) 
Qp{x) 



'4){x)ip{x) = L\x)R{x)+R'<{x)L{x), 
4){x)-t^ip(x) = R\x)R(x) + L\x)L{x), 
Vi(x)7V(x) = R\x)R{x) - L\x)L{x), 
tp{x)i-f5^p{x) = iL\x)R{x) - iR\x)L{x) 



(20) 

(21) 
(22) 
(23) 



6 



) become 



^s{x) = V^(a;)(7-i75V'(.7-) = -L{x)R{x) - R{x)L{x), (24) 

evix) = V^(x)C-i757X^) = -R{x)R{x) - L{x)L{x), (25) 

^^{x) = 4>'^{x)C-^^5l^Mx) = -R{x)R{x) + L{x)L{x), (26) 

^p{x) = '4F{x)C-^ii){x)=-iL{x)R{x)+iR(x)L(x), (27) 



and 



ls{x) = i'ixh^Ci'^ix) = -RHx)LHx) - L\x)R\x) ^ {is)\ (28) 

ty{x) s ^(x)7075C^^(x) = -L\x)L\x) - R\x)R\x) = (29) 

ty{x) = i>{x)i^^,Ci,^{x) = L\x)L\x) - R\x)R\x) = (ei.)^ (30) 

^p{x) = ^{x)iCtp^{x) = iR^{x)L^x) - iL^{x)R'<ix) = (Cp)^ (31) 

The reason of the identification of the symmetries of scalar S, vector V and pseudoscalar P of the Cooper pairs will 
be clarified later in this section. If the fields R and L are spinless, the vectorial Cooper pairs £_y and £,y are forbidden 
by the Pauli principle, gs and gp break the U{1) chiral symmetry and keep the U{1) gauge symmetry, while and 
are chiral invariant and break the U{1) gauge symmetry. Cj is given in a quadratic form of these densities in 
general. 

To understand the physical characters of ordered states, we examine the transformation properties under the charge- 
conjugation C, spatial inversion V and time reversal T symmetries. Especially, the broken T is interesting with respect 
to the relation with spontaneous magnetization of the system. The parity transformation in (l+l)-dimensions is given 

by 



tpixojXi) j^^'ip{xo,-xi), iIj{xo,xi) ^ ip(xo,-xi)-/°. (32) 



Therefore one finds 



esixo, xi) ^ gsixo, -xi), gyixQ.xi) ^ gyixo, -xi), 

gy{xo,xi) ^ -gy{xo,-xi), gp{xo, xi) ^ -gp{xo, -Xi), (33) 

^sixo,Xi) ^ ^sixo, -Xl), ^y{xo,Xi) ^ ^yixo, -Xi), 

^y{xo, Xl) ^ -^y{xo, -Xl), ^p{xo,xi) ^ -^p{xo, -Xi). (34) 
Hence, gs and ^5 are parity-even. Under the charge conjugation, 

Qs gs, Qv -Qv, Bv -Qv, Qp ^ -Qp, (35) 
^s^-b, C^^^"^ ei'^^'J., ^p^-^p. (36) 
The time-reversal transformation T is defined by 

V'(a;o, Xl) i 7°V'(-a;o, a;i), V'(xo, a;i) ^ tpi-xo, xi)7°, (37) 

with the rule to take the complex conjugation to matrix elements ( anti- unitary operations ). The densities will be 
transformed as 

gsixo,xi) ^ g*s{-xo,xi), gy{xo,xi) ^ gy{-xo,xi), 

gy{xo,xi) ^ -gy{-xo,xi), gp{xa,xi) ^ -g*p{-xa,xi), (38) 
^s{xo,xi) ^ Csi~xo,xi). ^y{x„.xi) ^ ^y{-xo,xi), 

^y{xo,Xi) ^ -^y{-XQ,Xi), ^p{xq,Xi) ^ -^p{-Xo,Xi). (39) 



In considerations on superconducting order parameters, one has to find both its transformation property under 

several groups and the Pauli principle [20]. In the case of this paper, especially we should consider the spin and the 
(1+1) -dimensional Lorentz ( namely, spinor ) symmetries. ( We must not confuse the "spin" and "spinor" degrees of 
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freedom. The spin degree is introduced as a kind of isospin in our model, while the spinor degree has the origin in 
(l+l)-Lorentz symmetry. ) The superconducting order parameter will be defined as follows: 

A(a;) ~ ^{x)i)'^{x), A{x) ~ i)'^{x)i){x), A = 'y°A'<'y°, A ^ CAC-\ A CAC'^ (40) 

Here, A and A are proportional to Gq; and we omit it for the notational simplicity of the discussion in this section. 
A and A have 4x4 matrix structures coming from the direct product of the two-dimensional spin and the two- 
dimensional spinor spaces. The (l+l)-dimensional Lorentz transformation is determined by the consideration of 
Lorentz algebra in the following form: 

S = expels), u = eoi, ugR, S-^=j°S^j°. (41) 

The definition of is crucial to keep the quantity like tijijj as a scalar. Under the (l-l-l)-dimensional Lorentz 
transformation, the superconducting order parameter will be transformed as 

A ~ VV-"^ ^ Sipij'^S'^ = e*T= VV''^e^T?' = i^ip'^ + | (V'V''^75 + TsV'V''^) + 0{w'^). (42) 

Hence, expanding A by the gamma matrices in (l+l)-dimensions as a complete set in 2 x 2 matrix spinor space, one 
finds 

A = (As + A^7'' + Apij5)c-^15, A = 75C(A5 + A^^t'' + Ap(-i75)). (43) 

Here, S ,V and P denote scalar, vector, and pseudoscalar, respectively. The scalar ASC75 has been chosen so as to 
be a (l-l-l)-dimensional parity-even order parameter under (32). ( The correspondences between the components of 
A and ^<j, £_y, are obvious for us: ^5 ^ AsC~^75, so forth. ) On the other hand, under the rotation of SU{2) 
spin space, the Dirac field is transformed as ( Tj ( i = 1, 2, 3 ); the Pauli matrices ) 

expQ6'iri)v, Oi G R, rf = -T2TiT2. (44) 
The order parameter matrix is transformed as 

A ~ VV'^ ^ e^^-^-VV'^e*^-^-^ = VV'^ + ^0i[Ti, ^Vt2]t2 + 0{e1). (45) 
Therefore, we find the following expansion appropriate to define spin singlet and triplet Cooper pairs: 

A = (A«+AfVi)zT2, A= -zT2(AW+AfVi), (i = 1,2,3). (46) 

Here, A*^^^ and A^'''' denote order parameters of spin singlet and triplet Cooper pairs, respectively. The spin singlet 
pairing of any Lorentz symmetry { S,V, oy: P ) becomes a scalar while the triplet has the vectorial nature under the 
spin rotation: Namely, the spin ( or, flavor ) and spinor ( Lorentz ) degrees of freedom are completely decoupled in 
our model. Due to the Pauli principle, the matrix A must change the sign under the transposition in spinor(g)spin 
space: 

A'^ = -A, A^ = -A. (47) 

By taking into account both the spin and Lorentz ( the (l-l-l)-dimensional spinor space ) degrees of freedom, we find 

the spin-singlet scalar pairing A^^''(a:)C75 ® ?T2 is allowed by the Pauli principle, and it corresponds to the s-wave 
BCS state. In this case A^ {x) = Ag {—x) ( an even function under parity ) has to be satisfied. On the contrary, 

the spin-triplet scalar ( a vector in spin space ) A^^''(a;) • C75 ® iTT2 has to become an odd-parity function under the 

spacial inversion to satisfy the Pauli principle: A^'{t) = —A^'{—x). If our model does not have the spin degree 
of freedom, only the pseudoscalar ( not the scalar ) component of the Cooper pair is allowed by the Pauli principle. 
A linear combination of ■ip'^C^^jcij''* ® iT2tJj and i'^^^'^^fi (E) iT2tp'^ corresponds to the Majorana-like mass term as 
RiT2R + R'^iT2R'^ and LiT2L + LHt2L\ and it cannot exist if our model is spinless. In (3-1-1 )-dimensional field theory, 
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a superconducting order parameter is a specific linear combination of ( parity breaking ) Majorana mass terms [23,24]. 
Usually, the symmetry of Cooper pairs will be defined in the momentum space, while A(a;) of a homogeneous system 
cannot have a momentum-dependence and hence it cannot become an odd-parity function in momentum space, thus 
spin-triplet Cooper pairs must be handled by the generalized BCS framework. 

Based on the discussion given above, we prepare the following interaction for generating both the s-wave dynamical 
Dirac mass ( coming from a nonzero VEV of gs ) and spin-singlet s-wave superconductivity ( again, from a nonzero 
VEV of ) from Cj in our theory: 

>Cf = ^i^i'f = 1^ W)' + ^i4'^r27,Ci,^){^^C-'J,iT2i>). (48) 

Here, we only keep the relevant terms for the dynamical generation of Dirac mass in a particle-hole channel and 
scalar ( in (l-l-l)-dimensional Lorentz space ) Cooper pairs in a particle-particle channel. The other terms appear 
in the Fierz transformation (??) may act as quantum fluctuations around a stationary point in the ( path integral 
) quantization of our theory, and they are suppressed by the large- A'^ limit because propagators of the fluctuating 
collective modes are proportional to 1/7V inside the effective action of the theory. This interaction can only handle 
s-wave pairings. For the generalization of p-, d- and /- wave pairings, especially to consider spin-triplet Cooper pairs, 
we will give the generalized BCS framework in Sec. IV. 

Another important aspect of our model under the context of the purpose of this paper is the Pauli-Giirsey ( PG ) 
symmetry, the symmetry of a rotation in particle- antiparticle space [51-53]. The transformation of the PG symmetry 
in (l+l)-dimensions is determined as follows [51]: 



PLi^ + Pr^ =(y f * ) , i^^i^PR + rPL = (i^ R), 



r ^ Pl^ + PrV = ( ^* ) . V'^ ^ i'^PR + i^PL = i-L, R^), (49) 



where. 



Pl = ^(1-75), Pr = ^(1 + 75), Pl+Pr = 1, PlPl=Pl, PrPr = Pr, PlPr = PrPl=0, 

= C^F y ^= -^^C-i = i-L, R). (50) 

Here, the definition of the charge-conjugation of ijj has been given in terms of R and L for the convenience of our 
discussion. ( More compact expression will be given in the form of displacement operators: 

Q = Q = {^,i,), Q' = exp[isF]Q, Q' = Qexp[-ieF], F = ^ , Q'Q' = QQ. (51) 

The quantity QQ, can be called as a "norm" or a "inner product", is "conserved" under the PG transformation. 
) Hence, in our definition of ip in this paper, the PG transformation corresponds to the additions of the complex 
conjugates R* or L* to a specific chirality sector i? or L of i/j. In order to examine the condition of the realization of 
the PG symmetry in our model, we rewrite the kinetic term Cq in the following form: 

-Co = -ipa (7°{«^o + vphp + m(^2)} + ivp^^di - m(fc2) j Va 

= {l°{ido - VFkp - H{k2)} + ivFl^di - m(fc2)) = ^y- (52) 

By performing the PG transformation defined in (49) to Co = {Cq + Cq)/2, we find that the Lagrangian Co has the 
PG symmetry at vphp + /^(fe) = 0, which is the same with the Lorentz symmetry condition for Cq, while the mass 
term m{k2) is neutral imdcr the PG transformation. The points as solutions of the equation = vpkp + /x(fc2) in the 
(fci, A;2)-surface of the energy eigenvalue may exist under the condition \vFkF/{2t)\ < 1: At = cos~^ {vrkp / {2t)), we 
cannot distinguish between tp and tp'^, 'ij) and tp'^: In this case, for example, R'^L and RL are the same composite field. 
We should note that this PG-symmetric condition has been obtained in the case of "free" field theory, and interactions 
( namely, renormalization ) between particles would modify it in general. The densities will be transformed under the 
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PG group as follows: 

(5(V;V) = ^Pr^' + i^PLi), S{'^j°^) = V^tV, ^CV-l V) = V'7 V, 5(V'75V') = ^Pr^'" - i^PLi>, 
<5(^75CV;^) = 0, (5(Vi7°75C^^^) = -V'^PflV' - 

^CV^tSsCV^^) = -V'^^iiV' + V^-PrV^^, Si^pCi;'^) = 2ipPRip. (53) 

Prom these results, we find that the four-fermion interaction {ip'^C~'^j5ip){'iljj5C'ip'^) is explicitly PG invariant ( 
ip'^C~^j5ip and ■ip'y^C-ip'^ arc PG singlet ). The displacements of the scalar S{ipip) and the pseudoscalar S{ipj5ip) show 
Cooper-pair-type fermion-fermion pairs. On the contrary, the vectorial fermion-antifermion ( fermion-fermion ) pairs 
give again fermion-antifermion ( fermion-fermion ) pairs under the PG transformation. 

Before closing this section, we should also mention the Vafa-Witten theorem [55]. This theorem states that the 
isospin ( flavor ) symmetry cannot be broken in Nambu— Jona-Lasinio-typc four-fcrmion-interaction models through 
the chiral condensation. However, this theorem was proved by using a set of Lorentz invariant eigenvalues by solving 
the two-flavor Dirac equations ipq = Xq { q = u,d quarks ). Because the Lorentz symmetry is explicitly broken 
in our model, the theorem could not be applied to our theory naively. For example, we consider the possibility for 
spontaneous generation of the vectorial condensate i^a{<^i)ai3ipi3 ( a kind of SDW ) which could break the rotational 
symmetry in the spin space of our model. Furthermore, the existence of the external magnetic field B will explicitly 
break the spin symmetry in our theory. 



IV. THE EFFECTIVE POTENTIAL AND THE GAP EQUATIONS 

In this section, we will obtain the eS'ective potential and the gap equations of several order parameters. By taking 
into account the results in the group theoretical consideration on the order parameters given in the previous section, 
we only consider the scalar components of the (l-l-l)-Lorentz symmetry. By using local ( a{x), Ag^\x) ) and bilocal 
( Xs{x,y), 'Es{x,y) ) auxiliary fields [56,57], the Lagrangian becomes 

iV 2N 



-YXs{x,y)V \x,y;x',y')xs{x',y')-^Es{x,y)W \x,y;x',y')Es{x',y') + ^^{x)M{x,y)'^{y). (54) 

Here, Ag' is the scalar, ,s-wave (iso)spin-singlet component of the order parameter A. Xs{x,y) and Es{x,y) are 
scalar, proportional to the bilinear forms V'(a;)(l -t- T)i'{y) and ^'^ {x)C~^-f5{l + T)iT2ip{y), respectively. They include 
both spin-singlet and triplet pairs in general, a and xs have no i7(l)-phase degree of freedom. In the expression of 
the Lagrangian, we have introduced the eight-component Nambu notation [16,22-25] in the following definition: 

*(^) = f &S V *(^) = (Mx), €{x)l (55) 



Vo^ix) 



and the matrix M. is defined as follows: 



where. 

Mil = l\ido + vpkF + m(^2) - ^^B^T-i} + ivpj^di -a- 771(^2), (57) 



M12 = A^s\x)C-^j5 «T2, M2I = -A'i\xy/5C iT2, (58) 

M22 = 7°^{ido - vpkp - l^{k2) + ^^BsTs} + ivpl'^di + a + m{k2), (59) 

Xsix,y) = xs\x^y) + xs\^^y) -^^ (60) 

Ss(a;,7/)= (E''g\x,y) + Ef\x,y)-TjC^^-/5(SiiT2, ^{x,y) = -iT2 (i$ -f^C (E'-g\x,y) + E'-g\x,y) ■ rj. (61) 
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Here, we have presented M in the form of Fourier-transformed expression in a;2-coordinate. x^s^ S^^' denote 

Is) (3) 

the spin singlet pairs, while Xs ^^'^ ■^s triplet pairs. The Fourier transforms of bilocal auxiliary fields and the 

nonlocal interactions become 

Xs{x,y)= [ /e-*^^-^«'-xs(P,g), Es{x,y) = [ /e-^^^-'«'^Ss(P,g), E^{x,y) = [ /e^^^+^'^'-S^(P,g), 

JpJq JpJq JpJq 

V{x,y;x',y')= [ [ / / e-'^^-^^'^'-^'"--^«''-V(P,g;P',<z'), 

Jp Jp' Jq Jq' 



W{x,y;x',y')= [ [ [ [ e-'^^-*^'^'-*«'^-*«''-V(P, P', g'), 

Jp Jp' Jq Jq' 



X ~^ y X -\- y I I I /re\\ 

= ^— , X= — - — , r = x-y, r=x-y. (62) 

Here, we have used the abbreviation defined in Eq. (13). X = {Xq,Xi,X2) denote the center of mass coordinates 
for fermion-antifermion and fermion-fermion pairs, while r = (xq — yo,xi — yi,X2 — yi) are the relative coordinates 

between particles they would form a pair condensation. Hereafter, wc take P^^ = P'^ = 0, choose the total momentum 
of any Cooper pairs at the origin of (/ci, /c2)-space. The partition function Z considered in this paper is defined as 
follows: 

Z = j I>^DV'2^crDA^^^r'A^^^r'XsI>Ss25Si'exp[i^ j dxodxiL . (63) 

It seems needless to say, however, we wish to emphasize that, this functional integral contains all of the information 
on the dynamics of the system. We employ the steepest descent approximation to the integration of the auxiliary 
fields in Z. The effective action F is found to be 

r = - ^ / dxodxi ) + 2 I dyodyidx'odx[dy'ody[[xsV~\s + ^sW-^Ei 



-UnBet^yM{x,y). (64) 

Here, Det^y is the determinant of both x and y spaces in the functional sense. We assume a and A^^' are 
constants, do not depend on spacetime coordinates. Hence the effective potential is obtained by its definition 
Veff = -r / {En J dxodxi): 

= ^ + + I jJ^{xs{k)V-\k,k')xs{k')+^s{k)W-\k,k')'Es{k')) + '-tr j^lndetM, (65) 

where the determinant becomes 

tr / IndetA^ = / ln[fco' _ (E+f] [k^ - {E^f] [k^ - {Elf] [k^ - {Ezf] . (66) 

Jk Jk 

The determinant has totally eight eigenvalues ±E^, ±El, ±i?+ and ±EZ if there is no degeneracy. Hereafter, we 
will take the static approximation, assuming that the order parameters xs, said the non-local interactions V and 
W have no ko- { energy ) dependence. This approximation corresponds to the neglection of retardation effects in the 
interactions. 

Due to the matrix structure of A^, it seems difficult to evaluate detA^ with including all types of order parameters. 
We will consider several situations: 

• Example (i): P3 ^ 0, cr ^ 0, x^s^ 0, Xs^ and A^g^ = E^g^ = = ( no superconductivity ). We assume 
Xs^ = (0,0,Xs3)- In this case, the energy eigenvalues become 

E+{k) - V'^fc? + i<^ + Mk2) + xi'\k) + xi(fe))2 T {vFkp + fi{k2) - ^B^) , (67) 

Eg (fe) = ^Jvj,kl + {a + m{k2) + xf {k) - xfi{k)Y T (vpkp + ti{k^) + ^^3) • (68) 

By this example, we could examine the magnetic-field induced chiral symmetry breaking ( chiral condensation, 
CDW/SDW-types ) in our model. 
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• Example (ii): ct 7^ 0, Xs'' + 0' ^ 7^ 0' ^ 7^ 0. ^ + 0' and B3 = Xs^ = 0. We will choose the spin-triplet 
Cooper pair as = (0,0, S^j). This choice is a special case of the condition for unitary order parameter 

of spin triplet Cooper pairs S^^S^^ oc 1, usually be employed in a (3+l)-dimensional theory [20] and in that 
case it indicates the time-reversal invariance of a spin-triplet superconducting state. The energy eigenvalues are 
obtained as follows: 



Etik) = J {yjvlkl + [m(fc2) +a + xf{k)? T [vfUf + n{k2)]f + |A« + + sg](fe)P, (69) 



Eiik) = y {^P^^ + [m(fc2) + a + x's\W T [vrkp + /.(fc^)])^ + |A^^^ + 4^^(fe) - ^'^l{k)\\ (70) 

EX and E~ are quasiparticles coming from positive energy states, while E'^ and E_ are quasiparticles coming 

from negative energy states. Duo to the coexistence of A^^^ Sg'' and S^^g, the (pseudo)spin degeneracies have 

been resolved, and totally four branches of the quasiparticle excitation spectra exist in this case. When S^^3^ = 
in this case, these spectra recover the ordinary BCS-type energy spectra: 



E^{k) = E+ik) = £±(fe) = y (^4fc2 + [m + a + xi'Hfe)]^ T [vrkp + /.])' + \A^^^ +E<i\k)\^. (71) 

Both E+ { = E+ = E^) and E_ { = Et = EZ) are doubly degenerate. Even if Xg ^ = A^^'+S^^^ = 0, the upper 
band E^{k) and the lower band —E_{k) have no intersection when m + a ^ Q. At ct -|- ■)^g\k) + 2ucos A;2 = 
with A^^^-hS^^^(fe) = 0, the energy dispersions become E± = Vp\k\\^ {vpkp + 2t cos ^2), and they correspond 

to ±5-1- discussed in the beginning of this paper. The pure-spin-triplet superconducting state 7^ 0, = 

= Ag^^ -|- 'E^p = 0, also gives similar spectra, have two doubly-degenerate excitation branches. 

We can construct the case where particles of up-spin component of tp participate the Cooper pair condensation 
while particles of down-spin component remain a Fermi gas [58] . In that case, the projection operator in the isospin 
space (1 ± T3)/2 should be inserted ( as vertices ) to particle interactions. 

Next, we employ the finite-temperature Matsubara formalism [59-61]. The /co-integration is converted into the 
Matsubara formalism by the following substitution in our theory: 

/^^X];^' ko^iojn, /3=^, (T; temperature), Un = ^—^^, (n = 0,±1,±2, • • •), (72) 

n 

and this method is convenient for our discussin in this paper. The four-fermion contact interaction model in (2-1-1)- 
dimensions is renornializablc [50], though the rcnormalization of the ofFoctivo potential could not be performed in 
the same way of the usual Lorentz-symmetric case. In fact, due to the existence of cosine terms in several functions 
of our theory, it seems difficult to employ the dimensional regularization, the standard method of regularization for 
usual (2-|-l)-dimensional field theory. Therefore, we do not consider the rcnormalization of the theory, and regard 
£ as a cutoff model. ( Integrals of (2-1-1 )-dimensional field theories usually become ultraviolet finite by employing 
the dimensional regularization. ) As our definition of the gamma matrices, our model is closer to (l-l-l)-dimensional 
fermion models than (2-1-1 )-dimensional field theories. Due to the Matsubara formalism, we should employ a non- 
covariant cutoff scheme. The effective potential Veff is obtained in the following form: 

(JTQ (-TO 

1 dki r dk2 dk[ r dk'^ 



''"2./_A 27r 2-K i_A 27r 27r 



{xs{k)V-\k, k')xs{k') + Es{k)W-\k, k')Es{k')) 

-\ £ ^ £ $ (^^ + ^+ + ^- + ^- + 1 + + + ^""""^^^ + ^"'"0 ' 

where, the cutoff A has been introduced to the fci-intcgral. The regularization of integration is controlled by A of 
fci-direction due to the asymmetry of the momentum space in our case. Because we use the linearized dispersion 
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approximation to fci-dircction to obtain our model from an underlying theory such as a Hubbard-type hopping model 
of particles on a two-dimensional lattice system, and we only take into account the long-wave-length components of 
the wave-fields of fermions move toward xi-direction, the relation A ^ u has to be satisfied. The cutoff A, the 
transfer t and u have mass dimension [mass]^, while vp is dimensionless. 

The stationary condition of K// under the variation with respect to order parameters ( i.e., the gap equations ) 
will be obtained as follows. The gap equation for a is found to be 







eff 



da 
2a 



1 dki r 
2y_A 27y_, 



dk2 



^tanhfi^i 



^tanhfi.; 



oa 2 



dE, 



da 



(74) 



while, for X5^(fe), 

dVeff 



= 



dxf{k) 



A 27r 



dk'n 



dEi 



dxi'\k) 



t&nh-E+ 
2 + 



dE+ 3 . dEZ 8 
tanh -EZ + ^= tanh -E~_ 

dxf{k) 2 dxf{k) 2 



^^tanh^E,- 
dxf{k) 2 



(75) 



The self-consistent equation for will be obtained by the replacement d/dx's' ~* ^/^Xss- On. the other hand, the 
gap equation for the spin-singlet s-wave Cooper pair A^^^ is found to be 

dVeff 



(1) 



,(3) 







2AW 
Go 



while, for Sg^^(fe), 



dki 



" dk2 
27r 



dE+ 



dE: 



p 



dA 



= tanh -El + tanh -EI + 

1) 2 + g^W 2 + 



dEt , p , 
: tanh -E+ 



dA 



(1) 



dEZ 



dA 



(1) 



:tanh-i;: 



(76) 



5S«(fe) 2 27r 27r ^ ^ > sK ) ^ 

dE+ p , dE- p 
-1-^^ tanh -EZ + tanh -E~_ 



dE+ p , 
tanh tiE+ 
9^(1) 2 + 



9^(1) 



^tanh^E,- 



(77) 



The gap equation for the self- consistent determination of E!3^^(A;) will be obtained by replacing d/d'B^-p d/d'E.^^g 
in the expression given above. Several derivatives appeared in the gap equations obtain different forms between the 
examples (i) and (ii). The derivatives appear in example (i) become 



dE^ 
~d^ 

dE^ 
~d^ 



dE^ 

dE^ 



dE^ 

dEl 



, (1) , (3) 

a + m + Xs +Xs3 



^vlkl + {a + m + xf +xfy 

, (1) , (3) 

a + m, + xs + Xsi 



dx 



(3) 



vlkl + {a + m + xf -X'si) 



m\2 



(78) 
(79) 
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The derivatives with respect to a and Xg ^ in example (ii) become 



dEX dEl 



da o^i^) 



vlkl + [a + xf + mf T [vpkF + iAT + I A^^) + 4'' + 41 ? V^^^i + + + xf? 



(80) 



dEg dEl 



^ vlkl + [g + x^'^ + mY T [vFkp + m] a + Xg ^ + m 



'{^v%k! + [a + x'i^ + mf T [vrkp + n]f + lA^ + _ 4) |2 ,J v%kf + {a + m + x^^^Y 
The derivatives with respect to ^ and S^g^ of the example (ii) become 

^-si ^ y (v/^fcf + + xi'^ + T [VFkp + /z])' + |A« + 4) + S^1|2 

_i A« + 4)(fc)-4)(fc) 



(81) 













dEl 


aA« 


d^ 

O^g 



^-si ^ J (^4fcf + [a + xi'^+m]2 T [vpkp + + | A^ + 4) - 4) p 



(82) 



(83) 



The dominant contribution lo [ lie iiitcgrjitioiis of the gap equations of the example (ii) will come from the quasiparticles 

of the vicinity of Vpkf + [a + X5 ^ + m(fc2)]^ — [vFkF + IJ'{k2)] ~ of the positive energy branch in the (fci, fc2)-space. 

Therefore, Vg// and the gap equations include Go, A, vp, u, t, gji-isB^ and T as external ( model ) parameters. 
We will calculate Vefj and the gap equations numerically in the space of these model parameters. In the case of 

kp = t = u = Xs'' = l^^s^ + ^ ± 1 = S3 = of the example (ii), the gap equation for a becomes that of the 
well-known result of the (l-l-l)-dimensional Gross-Neveu model by evaluating the method of large- A'' expansion [43]: 



^ = fj^E^f"''^'^^^^''^' E{k,)=E%{kuQ) = ^vlkl + a^ 



(84) 



Now, we specify the functional forms of V and W . Usually in the generalized BCS framework, the effective attractive 
interaction will be decomposed into channels by using the complete set of angular momentum eigenf unctions. In the 
quasi- (H-l)-dimensional case we consider here, the decompositions of the effective interactions V{k, k') and W(k, k') 
are simply performed by Fourier expansions in both fci and ^2 directions. Under the weak coupling approximation, 
V{k, k') can be written in the following form in general: 

V{k, k') = {gnm cos nki cos mk2 cos nk'^ cos mk'2 + cos nki sin mk2 cos nk'^ sin mk2 

n,Tn 

+9nm sill COS mk2 sin nk[ cos mk'2 + 9nm ''^^T- sin nk[ sin mk'^ . (85) 

Here, |fc| — \k'\ must be satisfied before and after scattering of particles. Parallel with this procedure, the order 
parameters xs and S5 also be decomposed as follows: 



Xs = IXo 



^ ( I Xnm I COS nki COS mk2 + \ Xnm I COS nki sin mk2 
+ \Xnm \ sinnfci cosmfe + \Xnm \ sinnfci sinmA;2y (86) 
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For example, we consider the following components of the partial-wave decompositions: 

V{k,k') = W{k,k') ~ go, (s), (87) 

~ sinfci sinfc^, (p^), (88) 

~ 5p„ sinfc2sinfc2, (py), (89) 

~ 5d^2_^2 cosfc2CosA;2, {d^2_y2), (90) 

~ 5d^^sinA;isinfc2sinA:isinA;2, {dxy), (91) 
~ Qf^ sin fci cos k2 sin fcj cos k'2, {fx)-, (92) 
~ ff/^sin2fc2sin2fc^, {fy). (93) 

Here, s, p^;, denotes the symmetry of fermion-(anti)fermion condensations, and we have identified that x = ki and 
y = k2- We should notice that, the parity of the order parameters will be determined by both the dependence of 
fci and ^2- Here, we assume V = W, while in general, we can choose different interactions between V and W. By 
putting these several types of interactions in (87)- (93), we will obtain the gap equations under the generalized BCS 

formahsm. 

We briefly comment on the renormalization of 14/ / . As mentioned above, due to the anisotropy of the model, the 
renormalization procedure combined with the dimensional regularization cannot be performed in our theory. Let us 
consider the case A = x = 2 = for our example. In the effective potential evaluated by the cutoff regularization, 
we can incorporate the fohowing renormahzation condition [43,62]: 



eff 



Gr 



(94) 



Here, Gren is the renormalized coupling constant, jXr denotes the renormalization point. Therefore, ^4// can be 
written formally as 

''^'^'^ o^ + m, /(.).K//-^. (95) 



Gren 2 Oa 



However, it is unclear for us that whether the cutoff A is removed or not by introducing the renormahzation point in 
the procedure given above. At least at the level of analytical formulation, it seems impossible. Moreover, we have to 
solve the renormalization group equation for obtaining the results they are independent on fir- When we consider the 
case cr 7^ 0, A^^"'^'' 0, xs ^ a-nd S5 ^ 0, the Coleman- Weinberg type renormalization condition will become more 
complicated one than Eq. (94) ( though we can write it formally ), and thus, in fact, it it impossible to employ the 
renormalization procedure for our 14//. 



V. SUMMARY 



In this paper we have constructed the BCS and generalized BCS theory for the dynamical mass and superconducting 
gap generations in the quasi-(l-|-l)-dimensional rclativistic model. We have introduced a "relativistic" ( Dirac ) 
model where particles move intrachain direction as a homogeneous gas while the motion of particles toward the 
interchain direction is described by a hopping term similar to the case of Hubbard model. We have performed the 
group-theoretical consideration on several fcrmion-(anti)fcrmion pair functions, and have chosen the suitable order 
parameters for the dinamical masses and superconducting gaps in the possible situation of condensed matters. We 
have constructed the effective potential and the gap equations in analytical manner. 

As the next work, we will present the results of numerical solutions of the gap equations. The thermodynamic 
properties, response to external fields will also be given. Preparation for the presentation of these results is now in 
progress, and we will publish them as the part II of this study. 
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